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Abstract 



Recursion relations for integrals of amplitudes over the phase space, i.e. for partial 
wave amplitudes, are introduced. In their simplest form these integrals are propor- 
tional to the s-wave amplitudes and represent rigorous lower bounds on the total cross 
sections. The connection with classical field equations in D dimensions is established. 
Previous results on multiparticle amplitudes are easily reproduced. 
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The high-multiplicity limit of processes involving scalar particles has been stud- 
ied in several contexts p], |3|, f|, [5|, ||, |7]. Ordinary perturbation theory for the 
amplitude A(oo — > \) predicts a factorial growth, n\, which is inconsistent with the 



unitarity limits. Higher-order corrections at threshold || ^, [L0| exhibit an even faster 
growth with n, and the whole perturbative expansion breaks down for sufficiently 
high multiplicities. 

On the other hand, the nullification phenomenon, i.e. the nullification of 



A(G — > \) amplitudes at threshold in certain theories || TT|, [12|, [13], including 
the Standard Model with specific relations among the different couplings [[14], |l5|j , 
might be the signal that, under certain conditions, the perturbative expansion is still 
consistent with the unitarity limits: delicate cancellations among different graphs 
contributing to the amplitude might soften the n\ growth below the limiting case 
of the a/ttT one, so that the cross sections still satisfy unitarity bounds. The case 
of sinh-Gordon interactions in two dimensions is a quite interesting example of such 
cancellations, where the nullification survives at any energy and/or at any order of 
perturbation theory 



In this paper we are introducing some new recursion relations involving the inte- 
grals of the amplitudes over the phase space, 

a(n,w) = J f[d D Pl 5( P 2 -m 2 )6(p°)A( _ e , . . . , \)^(E /> - ?) W 



where w 2 = P 2 . These integrals are proportional to the s-wave amplitudes. Further- 
more, using the well-known Schwarz inequality, we obtain 

a(n,w) ~ J ft d D Pl 5{p 2 - m 2 )6( Pl °)\A(^ ^ . . . , x )l 6 > '^pf (2) 

where a(n,w) is the total cross section and V n (w) is the phase-space volume. 

On the other hand, recursion relations for tree-order amplitudes are, at least 
formally, easily understood. For instance the amplitude *4(oo — > \) can be written 

as 

« _ %Xp+i ^({ v ) M }Aoo) 
-VV-V = ?5 \oo!((E^ J°°¥ ~ P) 



(3) 



where 



**n\ ^j) {■> \) — ^/oo) y / e, '" , ^/\ J KJ 

and V = { .i , ) = oo, . . . , | = oo,...,\), 5Z\) — \} represents a given 

V V 
permutation of the final state momenta. The sum over n\, . . . , n p has to be understood 
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under the constraint n\ + n 2 + ■ ■ • + n p = n. The coupling constants X p correspond 
to an expansion of the interaction potential 



When all the momenta are equal and on-shell, the kinematical configuration corre- 
sponds to the threshold case, and the recursion relation of Eq.(|3|) has been studied 
extensively using the well-known technique of generating functions ||. 

In order to combine the recursion relation of Eq. (|3|) with the integration over the 
phase space in Eq. ([!]), we introduce the Laplace transform of the integrated amplitude 
a(n, w) as follows: 

a(n,P) = J d D Pe- p - p a{n,w) ; (6) 
which taking into account the overall 5-function, this can be written simply as 

„ n 

3(n, 0) = jtl d D Pi 5( Pl 2 - m 2 )9( Pl )e-^A( e , . . . , a) . (7) 



i=i 



Defining 

d 

a(n,f3) = —in\(d 2 — m 2 )b(n, (3) , = 



d/3» v ' 

and using Eq.(|3|), we have 

(d 2 -m 2 )b(n,P) = J2^r- E 6(m,0)6(na, P) . . . b(n p ,(3) . (9) 

p=2 P' ni,n2,...,n p 

We can now proceed by using the method of generating functions. Defining that 

oo 

#;r) = ^e->,/J) (10) 

n=l 

we find that the generating function satisfies the classical field equations given by 

(d 2 - m 2 )(f) = V\(j>) . (11) 
The above equation can be seen as a generalization of the well-known result that 



the generating function of tree amplitudes satisfies the classical field equations JL7 



integrals of amplitudes over the phase space can be calculated as well using the 
classical solutions in D dimensions. On the other hand, since we seek only solutions 
that depend on the measure of the vector Eq.(|TTD becomes an ordinary (non- 
linear) differential equation 

which under certain circumstances can be solved exactly. 

However hand we have to be careful when we consider Eq. ( T2 ) in connection with 
Eq . (|10"D . Since we are dealing with a system of differential equations the answer might 
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not be unique. This fact is reflected in the dependence of the solutions of Eq.fJTJ]) 
on r, which will indroduce an 'infinite' arbitrariness. This can be illustrated in the 
following example. Let us choose b n (x) = ne~ nx , which is the Laplace transform 
of a 5-function, S(t — n). They satisfy the following system of ordinary differential 
equations 

d 2 

(— - ljb n (x) = 6 ^2 b ni (x)b n2 (x) (13) 



which can be seen as a one-dimensional 3 theory. Defining <ft = J2™=i e nT b n (x) 
Eq.(|l3"D becomes 

(&- 1 )^ 6 "' 2 (14) 

with the solution 

0(x;r) = isinh- 2 (^±|^) . (15) 

It is now obvious that the existence of an arbitrary function of r, Xo(r), does not 
allow to calculate the 'amplitudes' b n unless some other conditions are specified. For 
instance, if we require that the inverse Laplace transform of b n (x) vanish for t < n, 
we single out the initial solution. Indeed this requirement applies also in the case of 
D > 1, since all physical amplitudes are non-zero only for w > nm. 

Having the solution of Eq.(ff2j), we can find the integrated amplitude as follows 

a(n, w) = w^/ 2 ————^-- / dp I u ((3w) (3 D / 2 a(n, (3) (16) 



2 D / 2 - l Q D T{^) m Jc-ioo 
where u — D/2 — 1, 1^ is the modified Bessel function and 



f 1 if D = 2 

nD -\ 27TV/IW) ifD>3 ■ (17) 



Several results can be reproduced using the above formalism. First of all expand- 
ing around (3 = oo, which in w-space corresponds to an expansion around threshold, 
w = nm, we get (D = 4) 

oo 

b(n, (3) = e -npmp-a(*-W ^ b k (n)p~ k . (18) 

fc=0 

The first coefficient, bo(n), satisfies the equation 

n \ 

(n 2 - l)b (n) =E" E r b (niMn 2 ) . . . b (n p ) (19) 

p=2 P' ni,n 2 ,...,n p 

with the well-known result 0, |TTJ. Notice that the leading term in the expansion 
around w = nm, is proportional to n\{w—nm)~ !1 2~ /r(|(n— 1)) which, taking (w—nm) 
to be constant, is not unitarity- violating as we expected, since it corresponds to 
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the non-relativistic approximation. In contrast, the saddle-point approximation for 
the integral in Eq.fllDD should be used when w — > oo, but / = w/nm is constant. 
Assuming that in this limit, b(n, (3) ~ b (n) exp(— mp(/3)), the saddle-point equation 
is given by 

This leads to the following recursion relation for the leading coefficient 

/V - l)6d(n) = E E & o(^i)&oK) • • • b (n p ) (21) 

p=2 P' ni,n2,.:,n p 

where we have replaced n 2 by (n 2 — 1) in order to make the recursion relation consistent 
for all n > 1. It is easy to see that Eq.(|21|) is nothing but the recursion relation 
for the lower bound on the amplitude given in references [|] and that (p(/3) ~ 
log(i^i(/9m)//3), which is given by the initial condition that 

b{l^) = 27[mK f m) . (22) 

Another advantage of our method is that one can study massless theories as well, 
for which the threshold amplitude is singular. This is possible because, in several 
cases, the singularity of the propagator is cancelled by the phase-space integration. 
For instance, for the 4 in four space-time dimensions, it is easy to show that the 
result is indeed finite for the integrated amplitude and that Eq. @ is suitable for the 
study of multiparticle massless amplitudes. 

The amplitude A(E — > \) can be studied in exactly the same way. The resulting 
equation is given by 

n \ 

(d 2 + 2q ■ d)c{n, f3,y)^J2^T £ & y ^ ® ^ 

p=l P- m,n2y,np 

where y = (3 ■ q/ (3 and 

- in\(d 2 + 2q ■ d)c{n, (3,y) = J d L 'Pe^* a 2 (n, w,P-q) (24) 

with 

„ n 

a 2 (n, w,P-q) = JY[ d D Pl 6(p 2 - m 2 )6(p°)A £ (U; e , . . . , \)^(£ /> " V ) 

(25) 

where ^4 e (U; e , . . . , ,\) is the amplitude for the process 

V V v 

0(g') + 0(g) -> 0(pO + • • • + 0(p ft ) . 

Of course for the physical amplitude we have to take the limit q' 2 — > m 2 , or 2P • q — > 
w 2 , so that there is a one-to-one correspondence between the physical amplitudes and 
the poles of the function c(n, w,P ■ q) in {w 2 — 2P ■ q). 
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In order to study Eg. (|2"3"D , we expand around threshold, which in /3-space corre- 
sponds to an expansion around f3 = oo. The leading coefficient satisfies the following 
equation 

n \ 

n(n-2E)c (n,E) = J2^ 1 E c (m, E)b (n 2 ) . . .b (n p ) (26) 

p=l P- ni,n2,—,n p 

where 60 satisfies Eq. fllTf) . Equation (|26| ) exhibits the nullification properties, i.e. the 
existence of a finite number of poles in (n — 2E) for certain theories [11], [12 1 . The 
next-to-leading coefficient, c\(n,E), satisfies the equation (D = 2): 

n(n - 2E) Cl (n, E) + {n 2 - n - nE)c (n, E) + 2(1 - £ 2 )— c (n, £) 

n \ 

= E ( c iK £)&o(n 2 ) + pc (m, £)6i(n 2 )) . . . b (n p ) (27) 

p=l P- ni,n 2 ,...,n p 



where 61 (n) satisfies the equation 



A 



(n 2 - 1)6! (n) +n(n-l)6 (n) = E E &iK)&o(n 2 ) . . . &oK) . (28) 

p=2 P- ni,n2,—,n p 

As a check of these equations we consider the case of the sinh-Gordon potential 
in two space-time dimensions, which is given by (including the mass term) 

2 

77? 

U(<f>) = — (cosh(A0) - 1) . (29) 

The result of the recursion relation Eq. fl2~T|) , taking for simplicity m = 1, A = 1 so 
that A p+ i = for p even and A p+ i = 1 for p odd, is the absence of any pole other 
than E = 1 for a(n,E), as expected, since we know that the nullification for this 
model is exact for any energy and not only at threshold [16|. This is a nice test for 
the validity of our formalism. 

Let us now study the high-n behaviour of the function c(n, (3, y) in the limit 
n — > 00, w — > 00, /=w/nm=constant. This can be done by studying the inverse 
Laplace transform in Eq. (|24|) . The saddle-point equation in /3-space can be written 
as: 

d fl F(P,y)=P, (30) 

where F ~ — log(c). The solution is given by F = mp(/3), where (p satisfies Eq . (|20|) . 
Writing c = co(n, E) exp(— mp(/3)), and keeping the leading terms in n, we find that 

(f 2 n 2 - 2fnE)c (n, E) = J2^ E c o(n u E)b (n 2 ) . . . b (n p ) (31) 

p=l P- ni,n2y,n p 

where we consider E to be of order n. It is interesting to notice that if 60 satisfies 
Eq.(^TJ) then Co(n, E) still has a finite number of poles in (nf — 2E) for the theories 
with the nullification property. The only difference from the threshold case is a trivial 
redefinition of the couplings, A — > A// 2 , and the fact that the poles have been shifted 



5 



from E = n/2 to E = fn/2, as they should. It is also interesting to notice that 
this result does not depend on the dimensionality of the space-time, which enters 
only in subleading terms, and that it is indeed true for the case of the sinh-Gordon 
potential in two space-time dimensions. Although the question of nullification beyond 
threshold is of capital importance, a complete study of this problem in the context of 
4-dimensional theories is beyond the scope of this paper and will be given elsewhere 

El. 

We conclude by emphasizing that the formalism we propose in this paper leads 
to recursion relations for integrated (partial-wave) amplitudes, which can be used to 
reconstruct the total cross section. The simple integral of the amplitude over the 
phase space, which is proportional to the s-wave amplitude, is a rigorous lower bound 
on the total cross section. This formalism is quite general and has the very inter- 
esting property that the generating function of the Laplace-transformed integrated 
amplitudes, at tree order, satisfies the classical equation of motion in D dimensions. 
All previously known results can easily be reproduced, using appropriate approxima- 
tion methods. Furthermore, this formalism enables us to directly relate tree-order 
unitarity, in the limit of high multiplicities, with the classical field equations, and 
offers the unique possibility to study multiparticle production above threshold. 
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